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1 Introduction 



QED is commonly recognized to be a very successful quantum field theory. Never- 
theless, despite the excellent agreement found between the perturbative predictions 
of the theory and some experimental results, QED is still a poorly understood the- 
ory. The most unsatisfactory feature of QED is the so-called Landau pole problem. 
When a test charge is added to the vacuum of QED, pair creation around this charge 
totally screens it and the resulting charge observed at large distances is zero. That 
means that QED becomes 'trivial'. Formally we need a cutoff to make it nontrivial. 
Lattice (nonperturbative) QED is supposed to resolve this problem. 

There are different lattice formulations of QED (e.g., with compact or non- 
compact treatment of the gauge field degrees of freedom and with staggered or 
Wilson fermions, respectively) possibly not belonging to the same class of univer- 
sality. If so, on a physical ground one has to decide which version of QED is realized 
in nature. Old, i.e. experimentally well-established physical results should finally 
be reproduced (at least qualitatively) from lattice computations. 

If we consider QED as arising from a subgroup of a non-abelian (e.g., grand 
unified) gauge theory, then we are led to a compact treatment, a point of view we 
want to adopt throughout this work. Concerning the fermion degrees of freedom 
we are dealing with Wilson fermions which - to our knowledge - have not 

yet been thoroughly investigated in case of QED. 

In a theory with Wilson fermions the chiral symmetry is broken explicitly. Pre- 
sumably, it can be restored by fine-tuning the parameters in the continuum limit. 
But, for nonzero lattice spacing we can expect only a partial symmetry restoration 
at some value of the hopping parameter (j3) by cancellation of the Wilson 

mass term with the ordinary mass term in certain vertex functions at zero momen- 
tum 0, |J. If so, we can approach the continuum limit and the chiral symmetry 
restoration along the line k c (/3). It is another question, whether there the chiral 
symmetry becomes explicitly realized or spontaneously broken. 

It is worthwhile to note that QED can also serve as a 'toy' model of QCD for 
the study of, e.g., the chiral transition. Therefore, even the confinement phase of 
QED is of interest, since the fermionic observables should behave similar to the case 
of QCD. 

In this work we investigate the phase structure of 4c? lattice QED with special 
emphasis on the chiral limit. We consider the standard Wilson action (WA) and a 
modified one (MA) which is derived from WA by introducing suppression terms for 
lattice artifacts into the gauge action. We discuss different techniques to establish 
the chiral limit. On the one hand they are based on the convergence rate of the 
conjugate gradient algorithm to invert the fermion matrix, on the other hand they 
use observables sensitive to the occurence of zero modes. Some preliminary results 
have been published in 0. These techniques are applicable to lattice QCD studies 
with Wilson fermions as well. 

We concentrate on the study of gauge invariant observables as a first step, be- 
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cause for gauge variant quantities such as the photon and fermion correlators one 
has to take care of gauge fixing ambiguities @, 0, §]. 

The organization of the paper is as follows. In the second section we define both 
actions and the fermionic gauge invariant 'order parameters' which we calculated. 
The third section is devoted to the discussion of the matrix inversion. The results 
of our calculations for the quenched fermions and for the dynamical fermions are 
presented in the fourth and fifth sections, respectively. The last section is reserved 
for the conclusions and the outlook. 



2 Actions and observables 

The standard Wilson lattice action Swa(U, ip, ip) for 4d compact U(l) gauge theory 
(QED) is 

Swa = S g (U)+S f (U,^). (1) 
Ineq.(p Sq(U) is the plaquette (Wilson) action for the pure gauge U(l) theory 

S G (l7)=/3-X;(l-cos0p), (2) 
p 

where (3 = l/g% are , and U Xfl = exp(i9 x ^), 6 XIX G (—11,71] are the field variables 
defined on the links I = (x,fi) . Plaquette angles 6 P = 6 X .^ V are given by 6 X .^ U = 

The fermionic part of the action Sf{U,iP,iIj) is 



f=l x,y s ,s'=l 
M = 1-K-Q(U), 

Qxy = E [W • (1 - • U Xfi + Syp-n ■ (1 + 7^)-' • ul-w] , ( 3 ) 

where M. is Wilson's fermionic matrix, Nf is the number of flavours (we take 
Nf = 2 for the case of dynamical fermions), and k is the hopping parameter. 

In our calculations we also used the modified (compact) action Sma with Wilson 
fermions 

Sma = S wa (U)+5S g (U) , (4) 

where the additional term SSq is introduced to suppress some lattice artifacts, i.e., 
monopoles and negative plaquettes. The notion of (DeGrand-Toussaint) monopoles 
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can be introduced as follows || : the plaquette angle is splitted into two parts 

9 P = 9p + 2im P , 9 P e (-ir,ir] , n P = 0, ±1, ±2, (5) 

where 9p is associated with the (gauge-invariant) 'electromagnetic' flux through 
the plaquette and np is the number of Dirac strings passing through it. The 
monopole charge within an elementary 3d cube c is defined as the net number of 
Dirac strings passing through the surface of c : 

K c = ±- E Op = - E n P (6) 

Z7r PGde PGdc 

where the correct mutual orientation of the plaquettes has to be taken into account. 
The suppression term in eq. (|J) can then be written as 

SS G (U) = X K ■ E \K C \ + Xp ■ E(l - sign(cos^)) (7) 

c p 

with and Xp as 'chemical potentials' controlling the suppression rate of the 
corresponding artifacts. We have chosen (A# , Xp) = (oo, oo) in order to obtain 
maximal suppression of monopoles and negative plaquette values. Then for the 
quantized theory the suppression terms can easily be rewritten as Kronecker deltas 
in the measure of the functional integral. 

In the case of the pure gauge U(l) theory the complete suppression of mono- 
poles removes the phase transition at (3q ~ 1.0 [10, O, and a unique Coulomb 



phase at positive /3's appears [JTTJ] . For staggered fermions no chiral transition 
is observed as shown in The additional suppression of negative plaquettes 

drastically improves the overlap with the lowest state, at least, for the case of the 
photon ]TT] . 



It is known, that at least in the strong coupling region for the standard Wilson 
action (in the confinement phase) the ordinary mass term and the Wilson mass term 
cancel within the pseudoscalar mass mp at some k = K c (f3) , so that quadratic 
terms in the effective potential vanish for the pseudoscalar field , and n c = 0.25 at 

In the weak coupling range (Coulomb phase) perturbative calculations indicate 
that the mass of the fermion becomes equal to zero along the line k c (/3) (for the 
free field theory k c = 0.125) §. 

Besides the average plaquette ( P ) we calculated the following fermionic ob- 
servables 

(M = ^ ■ (T^M-^g , (8) 
<^75V>> = ^-(Tr^-M- 1 )^, (9) 
<n> = ^•(Tr(^- 1 75 ^- 1 75 )) G , (10) 
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where ( )g stands for averaging over gauge configurations which in presence of 
dynamical fermions includes the fermionic determinant. V = L 4 is the number 
of sites. The advantage of the order parameter in eq. (JT0|) - the 'pion norm' [[OJ 
- is that it appears to be a very sensitive observable in the 'critical' region. This 
can be understood by considering the spectral representation of the fermionic order 
parameters. Let f n = f n (s,x) be the eigenvectors of M. with eigenvalues A„ , 
and g n = g n (s,x) be the eigenvectors of 75.M with eigenvalues /x n : 

Mf n = A n • fn , lbMg n = fin' 9n ■ (H) 

Then one can easily obtain a spectral representation of the fermionic order param- 
eters : 



G 



< n >^?^' (12) 

Evidently, an eigenstate of M. with eigenvalue zero is also eigenstate of 75.M . So, 
the presence of configurations which belong to zero eigenvalues of M. also gives 
rise to poles in II (for more properties of the fermion matrix see, e.g. [Hp. 



It is commonly expected that close to k c the minimal eigenvalue |A m j n | 2 of 
M} M. is small and 
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i 1 -.) ■ < 13) 



Therefore, the appearance of small eigenvalues Aj ~ can be a good indicator for 
the chiral transition. For more discussion on this point for the Wilson action see, 



e.g., [15|, |16[ and references therein. 

In our calculations we have chosen periodic boundary conditions for fermions 
in all four directions. It is worthwhile to notice here that in the case of quenched 
fermions the choice between periodic and antiperiodic boundary conditions is irrel- 
evant. Indeed, the change of the sign of all, say, time-like links in the 'last' time 
slice does not change the gauge action, but this corresponds to the interchange of 
periodic and antiperiodic boundary conditions for fermions [[[7 . 



In order to extract information about the photon mass we have measured the 
photon plaquette-plaquette correlator at non-vanishing momenta and for several 



values of the coupling parameters (see [IS, [11 



3 Matrix inversion 

The inversion of the fermionic matrix M. = 1 — k ■ Q(U) defined in eq.(^|) is 
equivalent to solving the matrix equation 
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M-X = <p. (14) 

for any given 'input '-vector ip. 

One of the most popular methods is the conjugate gradient method (see, e.g., 
19| , [20[] and references therein). This is an iterative method to solve systems of 



linear equations D ■ X = (p, where D is a hermitian n x n matrix and ip is an input 
vector. The convergence of the conjugate gradient (eg) method is controlled by the 
condition number £ = A max /A m j n , where A max and \ m in are the maximal and 
minimal eigenvalues of D (D = M.^M. in this case). We observed for our standard 
cg-method that at a large enough number of iteration steps N cg the average residue 
(R) = (R)(N cg ) behaves as 



(R) = C ■ exp(-a • N cg ) for N cg > N , 

a = In^r (15) 

independently of the distribution of eigenvalues Aj provided n is large enough. This 
kind of functional dependence of a on £ was analytically obtained for some other 
gradient methods (e.g. Chebyshev methods) p0[ . To check the behaviour (|H)D we 
generated D with different (uniform, gaussian, double-peaked,. . .) distributions 
of eigenvalues and given X min and A max . The components of the input-vector <p 
were chosen randomly with gaussian distributions. 

In Fig.l we show the dependence of the average (normalized to unity at N cg = 1) 
residue (R) on N cg for a gaussian distribution of eigenvalues Aj of the matrix 
D in the interval [\ m in', A mQ x] as an illustration. Broken lines correspond to the 
behaviour exp(— a ■ N cg ) . One can see that after some steps N cg > N eq.(pT5|) 
works very well. The behaviour practically does not depend on the choice of the 
initial vector X : we checked this for a gaussian distribution of the components 
Xo(i) and for different 'ordered starts' (e.g., X (l) = 1, X (i) = at i > 1). The 
parameters Nq and C depend on the distribution of Aj. 

This observation gives the possibility to use the convergence of the cg-method 
to estimate k c . This estimation, however, is biased because of the dependence of 
the parameters N and C on the distribution of eigenvalues. 

Obviously, the fermion matrix inversion is the most time consuming step in the 
numerical work. Therefore, a maximal improvement of the convergence of the con- 
jugate gradient algorithm is highly important, in particular near k c . We have 
achieved considerable improvement by the use of preconditioning methods, e.g., 
polynomial preconditioning [^TJ plus 'even-odd decomposition' |22|, |23|. Further- 



more, the appropriate choice of the start vector X allows to speed up the procedure 
additionally: 



X n 



la 



1=1 



i+Y,K l -Q l -v, (is) 
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where Iq is some integer depending on k. For smaller K-values (far enough from 
k c ) the hopping parameter expansion (|16| ) with large enough l already represents a 
good approximation of the solution we search for. This approximation is even less 
time consuming than a combined application with the cg-method. 

4 Quenched fermions 

The quenched approximation of the theory can serve to provide very useful infor- 
mation about the behaviour of the fermionic order parameters and to develop the 
numerical tools. Its comparison with the full theory is instructive for the elucida- 
tion of the influence of the fermionic determinant. Moreover, it has the advantage 
to be less time consuming than the full theory. 

As it was already mentioned the phase transition at (3q ~ 1.0 disappears by 
suppressing lattice artifacts, and the quenched modified theory has only one phase 
for nonnegative /3's : the Coulomb phase. Qualitatively, the behaviour of the order 
parameters we measured in the modified theory MA at (3 > turned out to be the 
same as in the standard theory WA within the Coulomb phase {(3 > (3o). 

Let us now discuss the possibility to estimate K c {0) in the Coulomb phase based 
on the convergence of the cg-method. It follows from the discussion in the previous 
section that for the inversion of Ai^Ai the number of cg-iterations required for 
the convergence to some small R < e will behave as (N^ 1 ) ~ 1 — k/k c as long as 
N cg ^> N Q and k ~ n c . Therefore, the study of the dependence (A r ~ 1 )(K) for 
different lattice sizes L (at fixed coupling f3 ) can give a rather reliable estimate of 
k c . Fig. 2a represents the dependence of (N~ l ) on k at (3 = 1.1 (Coulomb phase) 
for quenched WA (qWA) for different V = L 4 up to 20 . By linear extrapolation 
{Ncg 1 ) ~ > one can obtain reasonable estimates of n c , which, as we shall see, are 
consistent with that obtained by other methods. The volume dependence becomes 
weaker with increasing: L . With increasing: matrix size n = 4L 4 also the variance 
of N cg with respect to different statistical configurations decreases. An absolutely 
similar dependence of (N^ 1 ) on k we observed for the modified theory MA at 
different values of the gauge coupling. [] 

The comparison of the two theories WA and MA at (3 = 1.1 (both in the 
Coulomb phase) shows that the difference in the convergence is rather small, and 
the estimated values of k c are close to each other. Therefore, the influence of lattice 
artifacts in WA is already small at this value of the coupling. With increasing f3 
the estimated value of k c {(3) shows a weak dependence on (3 . 

Figs.2b,c show the dependence of (N~ l ) on the inverse lattice size L" 1 at 
various values of k in the case of the quenched standard Wilson action at (3 — 1.1 
(b) and the quenched modified action (qMA) at (3 = 0.8 (c) respectively. For 
this dependence we again see no qualitative difference between Wilson action in the 

the confinement phase of WA the number of required cg-iterations is highly increased and 
one obtains even better estimates of k c than in the Coulomb phase. 
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Coulomb phase and the modified theory and it appears to be 



(N~ l ) L = (iV-% + 



C(«) 



(17) 



L 



Time histories of different observables (e.g., ipip , ip^ip , IT ) can give a signal of 
approaching k c . The appearance of sharp peaks in their time histories at distinct 
values of (3 and k testifies the appearance of small eigenvalues of the fermionic 
matrix Ai and, presumably, could be interpreted as approaching the chiral limit. 

Fig. 3 shows the time histories of the pion norm II at /3 — for the stan- 
dard Wilson gauge action (confinement phase) at different «'s. Strong coupling 
arguments predict a 'critical' value of k at k c = 0.25 . Indeed, one can see that 
below k c time histories are rather 'smooth' and do not show any evidence for 
peaks. The fluctuations become more pronounced with increasing k , and finally 
at k c ~ 0.25 sharp spikes appear. It is interesting to note that at k > k c these 
spikes do not disappear but instead become even stronger. With increasing (3 
(but at (3 < Po) time histories have a similar structure with the only difference 
that the spikes appear at somewhat smaller values of the Wilson coupling. Thus, 
e.g., K c (f3 = 0.8) < K c (f3 = 0) . The time histories of ipip and ■07s'0 look similar; 
here the spikes have a smaller order of magnitude in accordance with the spectral 
representations. 

On the contrary at (3 > (3 for Wilson action (i.e., in the Coulomb phase) 
the form of the time histories drastically changes (Fig. 4). We do not find peaks 
of comparable amplitude (~ 10 4 ) but, nevertheless, in a 'critical' region of k , 
the histories show up stronger fluctuations than at smaller or larger values of the 
Wilson coupling. Taking into account that for every configuration the pion norm 
II is the arithmetic average of 4 • L A terms corresponding to contributions of 
4-L 4 eigenvalues one can conclude that rather small eigenvalues occured. A similar 
picture as in the latter case also holds for the modified gauge action. 

The increase of fluctuations near k c can be best represented by the statistical 
variances of the observables discussed, especially the variance of II should be 
most sensitive. The value which we call (renormalized) variance o" 2 (II) is defined 
as follows 



where IT is the value of ■ J2 X y ^>Vy^~xyl^>-^~yxl^) for the i th configuration 
and iV is the number of measurements. ('Sp' denotes the trace with respect to 
Dirac-indices) . This enables us to locate k c ((3) through the maximum of a" 2 (II) . 

Figs.5a,b show the variance of the pion norm a 2 (II) in case of the quenched 
Wilson action qWA in the Coulomb phase at (3 =1.1 (a) and for the quenched 
modified action qMA at f3 — (b). In both cases there is a clear signal for the 
chiral limit or transition at some k c ~ 0.14 (a) and k c ~ 0.135 for the case (b). 



a 2 (II) 



V 1 

44 ' n 



N •. N 

£(n*-n) , n.i^ 



(18) 
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For increasing lattice size L the peak at these values of k c becomes more strong 
and narrow. 

Contrary to this the average pion norm (II) itself (see Fig.8a-crosses) was 
found to have only a very weak size dependence. Moreover, (II) shows for qWA 
in the Coulomb phase and qMA the same dependence on k and it is practically 
independent of (3 at fixed k < 0.125 . 

At (3 < Po in case of the Wilson action (confinement phase) the value cr 2 (IT) 
is well defined only below k c (/3) . As it can be seen from the time histories for the 
pion norm the region k > k c ((3) {(3 < fio) is 'infested' by the small eigenvalues 
of the fermionic matrix, so it is practically impossible to define the average pion 
norm (II) and other fermionic order parameters as well in this region even with 
high statistics ( 10.000 measurements in our case ). 

Based on our quenched data we draw the following phase diagrams in the (/3, k)- 
plane : as already mentioned the quenched WA has a confinement phase for (3 < (3q 
and a Coulomb phase ((3 > f3 ), and these phases are separated by a 1st or 2nd 
order transition line at (3q (P3L 59|). Moreover, we have found a 'critical' line 



k = k c {(3) which decreases from k c ((3 = 0) = 0.25 to k c ((3 = oo) = 0.125 . These 
lines subdivide the (/3, ft)-plane into four areas which in the dynamical fermion 
case are expected to become distinct phases. 

The quenched MA has only one Coulomb phase for positive values of (3. The 
behaviour of our observables corresponds to that in case of the Coulomb phase of 
qWA. We are also able to locate the line k = K c (/3) which extends over the whole 
/3-range ( k c {(3 = 0) ~ 0.13 and k c {(3 = oo) = 0.125 ), dividing the phase diagram 
into two possible phases in the case of full QED. 

At (3 = the compact WA and non-compact QED |25| agree because of the 



common compact coupling of gauge fields to fermions. Since the values of k c (/3) 
differ for the quenched WA and quenched MA in the extreme strong coupling limit 
((3 — > 0) our modified theory MA might have different chiral properties at strong 
coupling in comparison with the non-compact QED with Wilson fermions (see [J24|). 



5 Dynamical fermions 

For the dynamical fermion case (full QED) we investigated the phase structure and 
the chiral transition for lattice sizes 4 4 , 6 4 , 8 4 and 12 4 . 

We have simulated both actions (standard Wilson and modified) using the Hy- 



brid Monte Carlo algorithm |26| . Throughout our simulations we have used = 20 
molecular dynamics steps followed by an acceptance/rejectance decision at the end 
of each trajectory. The length of the trajectories was choosen for most of our runs 
to be N5t equal to unity, resulting into acceptance rates of about 95 % at the 
considered coupling values. At a few values of the couplings we have considered 
several stepsize values 6r on 4 4 lattices, where we estimated integrated autocor- 
relation times for the plaquette operator from runs of about 10 3 trajectories. As 
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far as we can tell, increasing the length of trajectories up to NSr w 3 does not 
considerably decrease autocorrelation times, while a further increase leads to small 
acceptance rates and to a deterioration of the achieved MC statistics. 

Concerning the simulations with the modified action (§) we have chosen a unique 
suppression term 5Sg(U) in the form 

5S G (U) = A-£(l+tanh(£cos0 P )) , (19) 
p 

since the Kronecker deltas in the functional integral measure cannot be realized in 
the molecular dynamics part of the Hybrid Monte Carlo procedure in a straight- 
forward way. With an appropriate choice of the parameters A and B (B < 0) 
the term (|T9"| ) leads to a suppression of negative plaquette values. The parameters 
in eq. ([TJ]) were chosen to be A = 10 and B = —10. We have monitored the 
distribution functions of local plaquette values : In Fig. 6 we display the distribu- 
tion function of the local plaquette operator at the same values of couplings with 
(crosses) and without (dots) suppression term as an example. We see, that this 
choice of parameters leads to an almost complete suppression of negative plaquette 
values. Moreover, we made sure that in a simulation with a cold start monopoles 
never occur. 

In the simulations with suppression term the length of the trajectories was 
choosen to be almost a factor of 10 smaller in order to achieve reasonable accep- 
tance rates of about 90 %. It appears as if the constraint on the action reduces the 
available phase space considerably and allows only for smaller moves in the Hybrid 
Monte Carlo algorithm. 

Now let us turn to the question how the phase diagrams of the quenched theories 
given at the end of the previous section are changed when we allow for dynamical 
fermions. 

For the Wilson action WA we expect to have at least one additional 'upper left' 
and one 'upper right' phase which for simplicity we will refer to as the '4rd' (upper 
left) and '3rd' (upper right) phases (see FiglOa). We do not consider the question 
of possible additional phases at larger K-values ||. 

At ft ~ 1.0 and k < K\ ~ 0.18 we observe the same phase transition as in 
the quenched theory. Thermal cycles and time histories for the average plaquette 
show a clear signal of metastability at ft's close to fto and at different k < K\ . 
As well as in the quenched case this transition is driven by monopole condensation 
in the confinement phase. However, there is no such signal for ipip , which remains 
approximately constant at fixed k and varying ft around ft as measured at 
k = 0.06, 0.08, 0.10 . Therefore, at the moment we cannot conclude that there 
is a first order phase transition in this part of the phase diagram in disagreement 
with [27j. Probably the metastable states are due to the monopole loops wrapping 



around the torus. 

At ft < fto thermal cycles with respect to varying k showed a typical hysteresis 
behaviour and a metastability for both : the plaquette and ipip . In Fig. 7a we 
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display the hysteresis loop for tpip at /3 — 0.8 . Furthermore, the abelian monopole 
density strongly decreases with rising k (Fig. 7b). In both of the figures the full 
(open) symbols correspond to the runs with decreasing (increasing) k . 

For the standard compact gauge action case we therefore expect a 'horizontal' 
first order transition line from (/?, K c ) = (Pi ^ /3q, K\ ~ 0.18) to some 
(0, K e ) = (p 2 £ 0.15, K2 ^ K±) , the latter point assumed to be the 'right-lower 
corner' of the 4th phase (cf. discussion below). This transition is also driven by 
monopole condensation in the confinement phase. Yet we do not know what kind 
of transition line the phase boundary between the 4th and the confinement phase 
is, i.e. the 'horizontal' line continued from (/?, k c ) = (fa, Kq) to the point 
(/3, k c ) = (0.,0.25) already known from strong coupling expansion. 

The study of the (gauge invariant) photon plaquette-plaquette correlator T 7 (r) 
can serve to distinguish between the different phases. In the confinement phase of 
WA r~(r) decays very quickly with increasing r , so that for the photon mass 
m 7 the conclusion m 7 > can be drawn. In the Coulomb phase the behaviour 
of the correlator is consistent with the expected zero mass behaviour. In the 3rd 
phase our data also indicate the existence of a massless photon, but still more effort 
is needed to draw the final conclusion about the photon mass in the two 'upper' 
phases. 

In the following we want to discuss the chiral limit or transition from the 
Coulomb phase to the higher-K phase (s). 

As in the case of the quenched theory the study of the time histories of different 
observables, in particular, II , can give a signal of approaching k c . However, in 
the case of the theory with dynamical fermions the fermionic determinant strongly 
suppresses the small eigenvalues, and there are no such sharp peaks in the time 
histories as in the case of the quenched theory. Nevertheless, the study of time 
histories shows that by approaching k c (/3) they become much more noisy than far 
away from it, testifying the appearance of small eigenvalues of the fermionic matrix 
M (cf. @). 

Again the statistical variance should be a suitable parameter to express the 
fluctuations of the observables. 

The behaviour of the pion norm (II) and its variance in the full QED is similar 
to the behaviour of these quantities in the quenched theory. 

Fig. 8a shows the dependence of (II) on k at (3 — 1.1 for the standard 
Wilson action. One can see that the pion norm has a pronounced maximum whose 
position and value practically does not depend on the lattice size at large enough 
L (L = 8 and 12 in our case). At smaller values of L there is still a finite volume 
dependence. For comparison we show the data points (crosses) for (II) at the same 
values of j3 and k , but for a quenched simulation on a 8 4 lattice - this illustrates 
the influence of the vacuum polarization effects caused by dynamical fermions. 

2 In our preliminary publication H we conjectured a tachyonic behaviour, i.e. m 2 < , in the 
3rd phase. 
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For the variance cr 2 (II) of the pion norm (Fig. 8b) the signal in the full QED 
case is not as strong as in the quenched theory (smaller by a factor of ~ 6), 
because of the effect of the fermionic determinant. Nevertheless, the maximum of 
the variance 0" 2 (IT) of the pion norm shows a steady growth with increasing L . 
(We expect the errorbars shown in the figure to be somewhat underestimated (at 
most ~ 15%) since we have seen strong autocorrelations in our data for cr 2 (IT) , 
especially around k c .) We also observe a shift of the maximum of cr 2 (II) 'to 
the left' the more L is increased. Therefore, the variance of the pion norm seems 
to exhibit the characteristic features of a finite size scaling behaviour of a singular 
point on finite volume systems. 

We are led to the conclusion, that a chiral phase transition in the thermodynamic 
sense really exists. We might speculate that the corresponding singularity appears 
to be rather 'weak', i.e., other operators and their fluctuations, like the fermionic 
condensates and the plaquette operator do not show a singular behaviour, as far as 
we can tell for the given lattice sizes. As the variance of the pion norm itself can be 
interpreted as a higher derivative operator we speculate that the singularity may 
be viewed higher order phase transition. 

Based on our observations we conclude that the variance of the pion norm cr 2 (II) 
is a convenient 'order parameter' for determining the position of the chiral transition 
in the Coulomb phase. On the other hand the average pion norm itself, its time 
histories and (N~ ) can provide estimates of k c . 

Turning to the investigation of the phase structure of the theory MA (eq.@) the 
situation drastically changes in comparison with the standard theory WA. Similarly 
to the quenched theory in the full QED with artifacts suppressed we see only one 
phase at k < k c {(3) , and only one unique phase just above the line n c {(3) . More- 
over, we do not observe metastable states when crossing the 'horizontal' line, and 
no hysteresis shows up in thermal cycles. The average plaquette and the fermionic 
condensate, respectively, turn out to behave smoothly with varying n at fixed f3 . 
This has been checked for f3 = and f3 = 1.1 . The corresponding curves look 
very similar. So, there is no sign for a first order transition. 

We have also explored the small /3-region of the phase diagram above the chiral 
phase transition k c (/3) in the standard Wilson action formulation. As it was already 
noted strong coupling arguments suggest for (3 = and for /t-values above k c (/3 = 
0) = 0.25 the existence of zero eigenvalues of the fermion matrix. It is a valid 
question how this singularity is approached in the dynamical fermion simulation for 
positive values of /3, and for k > k c (/3), i.e. above the chiral phase transition line. 
In particular one might ask whether this singularity of the theory extents to finite 
/3-values and whether a corresponding unphysical phase of the theory exists. 

We have simulated the approach to small /3-values at a fixed Wilson coupling 
value k = 0.3 > n c (P) on a 4 4 lattice for various /3-values in runs of 1000 measure- 
ments, i.e. 4000 trajectories. 

In this region of the phase diagram we encountered acceptance rate problems 
in the Hybrid Monte Carlo algorithm. For this reason we have decreased the step- 



12 



size 5t to 0.025 while keeping the trajectory length fixed at unity - the resulting 
acceptance rate turned out to be about 84 %. 

With decreasing (3 we observe a rapid increase of the number of the conjugate 
gradient steps N cg needed for the fermion matrix inversion (Fig. 9a). A linear 
extrapolation of the average inverse number (N' 1 ) — > of conjugate gradient steps 
(dotted line in Fig.9a) can be consistent with a singular point around (3 ~ 0.15 . 
In Fig. 9b we display the time evolution of the pion norm II at decreasing values 
of (3 . The more (3 is lowered the more sharp spikes with increasing amplitude 
are observed in the time evolution of the operator, similar to that in the quenched 
theory (compare with Fig. 3). The appearance of such spikes makes it practically 
impossible to calculate a statistically reliable average value. 

We therefore conclude that the point (3 ~ 0.15 at k — 0.3 is part of a singular 
line (3 = P(k), k > 0.25 which separates the 4th from the 3rd phase. 

The 4th phase could be expected to be a parity- violating phase (see 0]). In case 
of QCD the existence of a parity violating-phase for k > K c (/3) was numerically 
based on the behaviour of the fluctuation of the pseudoscalar density (ipjsip) • 
We studied the distribution of ^75^ at {(3 = ; k < k c ) and ((3 = 0; k > k c ) 
in order to find a signal for parity-violation in the quenched approximation of the 
standard Wilson action on a 4 4 lattice. The distribution width increases roughly by 
a factor of 6 while crossing k c from below. This is due to the appearance of small 
eigenvalues of M. for k > k c leading to large fluctuations of ip^^ip seen as spikes 
in the time history. The shape of the distribution does not change, especially we do 
not find any double-peak structure for k > k c . One has to repeat this analysis on 
a larger lattice to control the finite size effects. Preliminary our data do not show 
parity-violation for the considered (f3, ft)-values. 

The phase diagrams for both theories (with standard gauge action and with 
modified action) are shown in Figs.l0a,b. For the standard action WA (Fig. 10a) 
we show four different phases at k^0.3 , whereas for the modified theory MA only 
two phases survive in the same K-range (FiglOb). 



6 Conclusions and Outlook 

Now let us summarize our results. 

• We have studied the phase structure of two theories with Wilson fermions and 
compact action with U(l) symmetry : standard Wilson theory and a modified 
one with lattice artifacts suppressed. Phase diagrams of both theories were 
shown. 

• For the standard Wilson theory there is a 'horizontal' line k c (P) from a 
(Pi, Ki) to (00; |) which separates the Coulomb phase from the 3rd phase and 
which corresponds to the chiral transition in this theory. Presumably, a higher 
order phase transition occurs on this line. At smaller values of P (P < Pi) 
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on the continuation of this line we find a presumably first order chiral phase 
transition. 

• When suppressing the lattice artifacts the latter 1st order phase transition 
disappears. 

• Nevertheless, after suppressing the lattice artifacts a chiral transition (i.e., 
appearance of near-to-zero eigenvalues of M. ) remains for all j3 values we 
considered. The statistical variance of the pion norm gives us a pronounced 
signal for both (quenched and dynamical) cases. Its location turned out to be 
in reasonable agreement with that obtained from the cg-method (measuring 
(N' 1 )). 

• The behaviour of all considered gauge-invariant observables along the critical 
line k c (/3) in the Coulomb phase is smooth. Thus, for the modified action, 
we have no indication for a (tri— ) critical point on this line, in contrast to 
lattice QED with staggered fermions and the non-compact gauge action. 

• The calculation of the (gauge- variant) fermionic propagator needs a gauge- 
fixing procedure which can appear to be rather nontrivial because of gauge 
fixing ambiguities 0. In a recent paper j7|, some of the gauge copies were 
shown to produce a photon propagator with a decay behaviour inconsistent 
with the expected zero mass behaviour in the Coulomb phase. In Ref. || it 
was proven that the "bad" gauge copies are due to the existence of pairs of 
Dirac sheets as gauge artifacts. For matter fields to be taken into account this 
problem requires additional studies before one can reliably compute fermionic 
correlators. 

• As a further step we plan to investigate the spectra of both theories and to 
determine renormalized quantities such as the renormalized coupling. 
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Figure captions. 



Fig.l The average residue (R) as a function of N cg for gaussian distribution of 
eigenvalues Aj of the hermitian matrix D in the interval [\ m i n ] X max ]. Broken lines 
correspond to the behaviour exp(— a ■ N cg ) where a is defined in eq.flLlf). 

Fig. 2a The dependence of (N~ l ) on k at /3 — 1.1 for the quenched Wilson 
action (qWA) and different lattice sizes L. The broken line has been added to guide 
the eye. 

Figs. 2b, c (iV" 1 ) as a function of the inverse lattice size L~ x at different values 
of k for (b) qWA at (3 = 1.1 and (c) the quenched modified action (qMA) at 
(3 = 0.8 . Broken lines are to guide the eye. 

Fig. 3 Time histories of II at different k at (3 = for qWA. 

Fig. 4 Time histories of II at different k at (3 = 1.1 for qWA. 

Figs.5a,b cr 2 (II) as a function of k for qWA at (3 = 1.1 (a) and the qMA at 
(3 = (b). The broken lines indicate the maxima of a 2 (II) . 

Fig. 6 Distribution functions of the plaquette operator without (dots) and with 
(crosses) suppression term fll9"|) at (3 = 0.8 , k = 0.1 on a 6 4 lattice. 

Figs.7a,b The thermal cycle for ipip (a) and for the monopole density p mon . (b) 
with respect to k at (3 = 0.8 for WA with dynamical fermions on a 6 4 lattice. 

Figs.8a,b Average pion norm (II) (a) and er 2 (II) (b) as a function of k for WA 
at (3 = 1.1 . Lattice size is 4 4 (circles), 6 4 (squares) , 8 4 (diamonds) and 12 4 
(triangles). For comparison the data from a quenched simulation on a 8 4 lattice 
is shown (crosses) in Fig.8a. 

Figs.9a,b (iV" 1 ) as a function of (3 at K = 0.3 (a) and time histories of II at 
k = 0.3 and and different (3 (b) for WA on a 4 4 lattice. 

Figs.lOa,b Phase diagram in the ((3, ft)-plane for dynamical fermions for the 
standard Wilson action WA (a) and for the modified action MA (b). 
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